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IMAGES IN A DIVERGENT SCANNING COMPUTED TOMOGRAPHY SYSTEM 

RELATED APPLICATION 

This appUcation claims the benefit of U.S. Provisional Application No. 
5 60/405,096, filed August 21, 2002, the entire teachings of which are incorporated herein 
by reference. 

BACKGROUND OF THE INVENTION 

The present invention relates generally to 2D and 3D computerized tomography 
(CT). In particular this invention relates to methods and systems for reconstructing 

1 0 proj ection data which are neither equlinear or equiangular in nature. 

In conventional computerized tomography for both medical and industrial 
applications, an x-ray fan beam and an equilinear or equiangular array detector are 
employed. Two-dimensional (2D) axial imaging is achieved. While the data set is 
complete and image quality is correspondingly high, only a single slice of an object is 

15 imaged at a time. When a 3D image is acquired, a "stack of slices" approach is 

employed. Acquiring a 3D data set one slice at a time is inherently slow. Moreover, in 
medical applications, motion artifacts occur because adjacent slices are not imaged 
simultaneously. Also, dose utilization is less than optimal, because the distance 
between slices is typically less than the x-ray collimator aperture, resulting in double 

20 exposure to many parts of the body. 



In a system employing true cone-beam geometry, a cone-beam x-ray source and 
a flat 2D equilinear or curved 2D equiangular area detector are employed. An object is 
scaimed, preferably over a 360-degree range, either by moving the x-ray source in a 
scanning circle around the object while keeping the 2D area detector fixed with 
reference to the source, or by rotating the object while the source and detector remain 
stationary. Jn either case, it is the relative movement between the source and object 
which affects scanning. Compared to the 2D "stack of slices" approach for 3D imaging, 
the cone-beam geometry has the potential to achieve rapid 3D imaging of both medical 
and industrial objects, with improved dose utilization. 

The cone-beam geometry for 3D imaging has been discussed extensively in the 
literature, as represented by the following: M. Schlindwein, "Interactive Three- 
Dimensional Reconstruction fi-om Twin Cone-Beam Projections", IEEE Trans 
NucLSci., Vol. NS-25, No. 5, pp. 1135-1143 (October 1978); Gerald N. Minerbo, 
"Convolutional Reconstruction from Cone-Beam Projection Data", IEEE Trans. 
NucLSci., Vol. NS-26, No. 2, pp. 2682-2684 (April 1979); Heang K. Tuy, "An 
Inversion Formula for Cone-Beam Reconstruction", SIAM J. Math, Vol. 43, No. 3, pp. 
546-552 (June 1983); LA. Feldkamp, L.C. Davis, and J.W. Kress, "Practical Cone- 
Beam Algorithm", J. Opt. Soc. Am. A., Vol. 1, No. 6, pp. 612-619, (June 1984); Bruce 
D. Smith, "Image Reconstruction from Cone-Beam Projections: Necessary and 
Sufficient Conditions and Reconstruction Methods", IEEE Trans. Med. Imag., Vol. MI- 
44, pp. 14-24 (March 1985); and Hui Hu, Robert A. Kruger, and Grant T. GuUberg, 
"Quantitative Cone-Beam Construction", SPIE Medical Imaging HI: Image Processing, 
Vol. 1092, pp. 492-501 (1989). 

Several methods for collecting cone beam data have been developed. One 
technique involves acquiring volumetric image data using a flat panel matrix image 
receptor, as described in U.S. Pat. No. 6,041,097 to Roos, et al. Another method uses 
image intensifier-based fluoroscopic cameras mounted on a CT-gantry type fi-ame. Such 
a system is described in a paper presented at SPIE Medical Imaging Conference on Feb. 
24 1997, by R. Ning, X. Wang, and D.L. Conover of Univ. of Rochester Medical 
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Center. 

U.S. Pat. No. 5,319,693 to Eberhard, etal. discusses simulating a relatively large 
area detector using a relatively small area detector by either moving the actual area 
detector relative to the source, or moving the object relative to the detector. 
5 However, there is a significant limitation of cone-beam reconstruction when 

individual flat detectors are reconstructed independently. Simply combining separate 
reconstructed portions of the object fi-om independently processed projections results in 
an image characterized by discontinuous jumps between the various projections. 
Alternatively, one could first combine the discreet data sets from each detector into a 
10 new single data set that is then reconstructed. However, by simply combining the data 
into a larger data array and performing standard reconstruction techniques, the data 
elements in the new data set are not equally spaced. Thus, the resultant images will be 
distorted geometrically, or the dynamic range of the reconstructed data set will not 
represent the true transmission values of the object being imaged. 

15 SUMMARY OF THE INVENTION 

The deficiencies in existing methods for combining image data from multiple 
flat panel detector arrays result from the fact that these detector arrays have neither 
equilinear nor equiangular geometries. The present invention relates to improved 
systems and methods for reconstructing projection data, including x-ray projection data 

20 for two-dimensional (2D) fan-beam and three-dimensional (3D) cone beam CT imaging, 
in which the geometry of the detectors is neither equilinear or equiangular, by 
reprojecting the actual measured data into a new virtual data array, which has an 
equilinear or equiangular geometry. In one aspect, multiple discreet projection data 
sets, which, when combined, are neither equilinear or equiangular, are reprojected into a 

25 new virtual data set on an equilinear spaced detector on a line or plane, or an 

equiangular spaced detector array on an arc or cylinder. The resulting virtual projection 
data set can then be reconstructed using standard backprojection techniques and 
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generate images which are geometrically correct, and represent the true x-ray 
transmission properties of the object being imaged. 

In one embodiment, the projection data from two or more ID linear or 2D flat 
detector arrays are reprojected onto a single equilinear or equiangular virtual detector 
5 array prior to filtering and backprojecting the projection data. 

In another embodiment, the projection data from two or more discrete detector 
positions are reprojected onto a virtual detector array having an equilinear or 
equilangular configuration, and the reprojected data is reconstructed to provide an 
image. 

10 The "virtual" detector array of the present invention is a data array comprising a 

plurality of pixels, having an equiliner or equiangular geometry, where the data values 
assigned to each pixel in the virtual array is based upon data from an actual detector or 
set of detectors having a non-equilinear and non-equiangular geometry. 

The present invention advantageously allows for the 2D and 3D tomographic 

1 5 reconstruction of objects. This invention enables divergent x-ray 2D fan beam or 3D 
cone beam tomographic reconstruction using a discrete number of ID linear or 2D flat 
detectors angled relative to one another by using a novel rebinning and reprojection 
technique onto virtual equilinear or equiangular detector arrays prior to performing 
standard filtered backprojection tomographic reconstruction techniques. 

20 The present invention is particularly usefiil for medical imaging applications, as 

well as numerous industrial applications, such as testing and analysis of materials, 
inspection of containers, and imaging of large objects. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The foregoing and other objects, features and advantages of the invention will be 
25 apparent from the following more particular description of preferred embodiments of 
the invention, as illustrated in the accompanying drawings in which like reference 
characters refer to the same parts throughout the different views. The drawings are not 
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necessarily to scale, emphasis instead being placed upon illustrating the principles of the 
invention. 

Fig. 1 shows standard equilinear and equiangular geometries used in various 
generations of CT scaimers; 
5 Fig. 2 shows a standard equilinear detector geometry in which the detectors are 

arranged with constant spacing along a line or plane; 

Fig. 3 shows the radiation profile of an imaged object defined by an equilinear 
arrangement of detectors; 

Fig. 4 shows a standard equiangular detector geometry in which the detectors are 
10 arranged with constant angular spacing along an arc or cylindrical surface; 

Fig. 5 shows the radiation profile of an imaged object defined by an equiangular 
arrangement of detectors; 

Fig. 6 shows three equilinear-spaced detector arrays positioned and angled 
relative to one another, resulting in a geometry that is neither equilinear nor equiangular; 
1 5 Fig. 7 shows the predicted radius of reconstructed object with three detector 

arrays generally positioned along an arc having a radius centered at the x-ray focal spot; 

Fig. 8 shows the predicted radius of reconstructed object with three ID linear or 
2D flat plate detector arrays positioned along a straight line in an equilinear 
arrangement; 

20 Fig. 9 is a flow chart diagram of the rebinning algorithm for reconstructing ID 

fan beam or 2D cone beam projection data which is neither equilinear or equiangular; 

Fig. 10 shows the projection of multiple angled detector array positions onto a 
single virtual flat equilinear detector array; and 

Fig. 1 1 shows the projection of multiple angled detector array positions onto a 
25 single virtual curved equiangular detector array. 

DETAILED DESCRff TION OF THE INVENTION 

A description of preferred embodiments of the invention follows. 
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Referring to Fig. 1, equilinear and equiangular detector geometries are depicted. 
A radiation source 13 projects radiation onto multiple one-dimensional linear or two- 
dimensional planar area detector arrays 14 that are angled generally along line or a 
circle. The detector arrays generate scan data from a plurality of diverging radiation 
5 beams, i.e., a fan beam or cone beam. The source and detectors are rotated around the 
object to be imaged, and a plurality of projection images is captured to computer 
memory for tomographic projection image processing. 

In the case of an equilinear geometry, a single source produces a fan or cone 
beam which is read by a Hnear ID or 2D array of detectors, as shown on the left. In the 
10 case of an equiangular geometry, such as shown on the right, the detectors occupy a ID 
arc to image fan beam data, or a 2D cylindriacal surface to image cone beam data. 

Referring to Figs. 2-3, an equilinear detector geometry is more clearly defined. 
As shown in Fig. 2, the detector elements in an array are arranged with constant spacing 
along a straight line or a flat plane. The angle between rays connecting the x-ray source 
15 point and the detector elements does not remain constant. A radiation absorption 

profile, or image, is generated with varying amplitudes for a region between the bank of 
detectors and the x-ray source, as shown in Fig. 3. Each ray is identified by its distance, 
s, from the projection of the central ray (s=0), and the absorption profile is denoted by 
the fiinction Rp(s). 

20 Figs. 4-5 illustrate an equiangular detector geometry. As illustrated in Fig. 4, the 

detector elements in an equiangular array are arranged with constant angular spacing 
along a circle or cylinder. In an equiangular geometry, in contrast to equilinear 
geometry, the angle between rays connecting the x-ray source point and the detector 
elements remains constant, but the distance between detectors may change. Fig. 5 

25 shows the radiation absorption profile for the region between the bank of detectors and 
the x-ray source. Each ray is identified by its angle, y, from the central ray, and the 
absorption profile is denoted by the fiinction RpCy). 

In many radiation imaging applications, it is desirable to image objects that are 
wider than the field-of-view of the detector array. One method for achieving a wide 
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field-of-view is to use multiple ID or 2D detectors, arranged end-to-end and angled 
relative to one another, as shown in Fig. 6. Another technique is to use a single array, 
translated to discrete positions along an arc opposite the x-ray source, to obtain a large 
"effective" field of view. In either case, when one or more equilinear ID linear fan 
5 beam or 2D planar cone beam detector arrays are positioned and angled along an arc 
opposite the x-ray source, the resulting geometry is neither equilinear or equiangular. 
As illustrated in Fig. 6, the projection of equally spaced detector elements, d^, on the 
angled arrays do not project onto equally spaced detectors, Pj, located on lines, planes, 
or arcs. For example, assuming the detector elements on the tilted arrays are equally 

10 spaced, the process of reprojecting these elements onto a new virtual detector array 
which is coincident or parallel to the central detector array will resuh in projections that 
are not equidistant. Hence, assuming the Fourier transform filtering is performed 
continuously along the axes of the angled arrays without resampling, the spacing of 
detector arrays cannot be assumed to be equal. 

15 Fig. 7 shows in more detail the resuh of filtering on non-equally spaced 

detectors. The predicted radius of a reconstructed object, R,, is calculated on an angled 
detector geometry, assuming fihering is performed without resampling onto an 
equilinear array. If we assume that the scanner focal length, FL, is 1000 mm, and the 
length of the detectors, L, is 400 mm, then 8, = 2*arctan((L/2)/FL) = 0.395 rad = 22.632 

20 deg., and R, = (FL/2)*sin 0^ = 192.31 mm. 

Fig. 8 illustrates this same calculation of the predicted radius of the 
reconstructed object assuming the same input parameters of focal length and detector 
length, but where the detector arrays are arranged in a plane to provide equilinear 
geometry. Here, 0, = arctan(L/FL) = 0.3804 rad = 21.795 deg., and = (FL/2)*sin 0, = 

25 185 .695 mm. The predicted radius of the angled detector geometry is larger than that of 
the equilinear detector geometry, and the resultant images with the angled detector 
geometry will be distorted geometrically. 

This problem can be overcome by reprojecting and resampling the data fi-om the 
angled detector arrays onto a "virtual" equilinear or equiangular array. The algorithm 



shown in Fig. 9 describes a method of reconstructing fan beam or cone beam x-ray 
projection data of an object, where the detector configuration is neither equilinear or 
equiangular. In particular, the algorithm describes a method for generating a new virtual 
equilinear or equiangular fan beam or cone beam detector array which is defined along a 
straight line or generally along an arc. For every pixel defined in the virtual detector 
array, the projection point in the original projection data is determined and the x-ray 
absorption amplitude for that point is calculated by interpolating the nearest neighbor 
pixels. Once resampling is completed, standard filtered backprojection and algebraic 
reconstruction techniques may be performed to generate image data. 
The method consists of creating a single virtual detector array for each projection 
position, which is defined as being equilinear or equiangular, and reprojecting two more 
real detector arrays onto the virtual array. Once the real projection data is reprojected 
onto the virtual detector, the data is filtered and backprojected using standard 
tomographic reconstruction techniques; 

As shown in step 101 of Fig. 9, the projection angle index, iproj, is first assigned 
the value 1 . At step 102, the x-ray source and detector array(s) are moved to a 
projection angle relative to the object being imaged. This can be accomplished by either 
moving the source and detector relative to a stationary object (preferably by moving the 
source and detector in a circle or arc around the object), or by keeping the source and 
detector stationary and rotating the object to the desired projection angle. 

The projection data can obtained for a plurality of projection angles (1 . . .nproj), 
preferably at a plurality of equally spaced angles as the source/detector and object are 
rotated 360 degrees with respect to each other. 

At step 103, a new virtual equilinear or equiangular array, P, is allocated. The 
virtual array, P, includes virtual pixels which are equally spaced in distance along a line 
or plane in the case of a virtual equilinear array, or equally spaced in angle along an arc 
or curved plane in the case of a virtual equiangular array. 

At step 104, the real projection data, D, fi-om each real detector array (1 . . .ndet) 
is acquired for the given projection angle, iproj. 
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For each real detector array, D, the real projection data is then reprojected onto 
the virtual array, P, at step 107. 

As shown at steps 108-1 15, the reprojection subroutine includes looping through 
each virtual pixel in the virtual array, P, (step 109), and for each virtual pixel, 

5 determining the real detector pixel, d, that is intersected by the line connecting the 
virtual pixel and the x-ray source (step 111). 

Once this actual pixel, d, is determined, an interpolation technique then is 
applied to d and its nearest neighbors on the real detector array to compute an x-ray 
absorption amplitude value to be assigned to the virtual pixel, p (step 1 12). This 

10 process is repeated until absorption amplitude values have been assigned to each of the 
virtual pixels in the virtual array. 

Once each of the real detector arrays has been projected onto a virtual equilinear 
or equiangular array, data from the virtual detector array is then filtered at step 1 1 7 and 
backprojected at step 118. As the name implies, there are two steps to the filtered 

1 5 backprojection algorithm: the filtering step, which can be visualized as a simple 
weighting of each Fourier transformed projection in the frequency domain, and the 
backprojection step, which can be seen as the dual, or in a more strict mathematical 
sense, the adjoint, of projection. Instead of projecting density values to a projection 
value, a projection value is backprojected, or smeared out, over the image points along 

20 the ray. This entire process is then repeated for each of the projection angles. 

Referring to Figs. 10 and 1 1, the process of reprojecting x-rays onto virtual 
equilinear and equiangular detector arrays is schematically illustrated. In Fig. 10, once 
actual projection images are captured, a new equilinear virtual detector is allocated and 
defined along a one-dimensional line in the case that a fan beam geometry, or along a 

25 two-dimensional flat plane in the case of a cone beam geometry. In Fig. 1 1 , the images 
are captured by the actual three-panel detector array, and then a new equiangular virtual 
detector is allocated. The equiangular virtual detector is an arc in the case of a fan beam 
geometry, and a curved cylindrical surface in the case of a cone beam geometry. In all 
of these embodiments, the new virtual array assiunes that detector elements are equally 
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spaced in distance or angle, respectively. For each detector element in the virtual array, 
the projected position in the real detector arrays is computed and an interpolation 
technique is applied to nearest neighbors on the real array to compute the correct x-ray 
absorption amphtude of the object to be reconstructed. Once the real detector arrays 
5 have been projected onto the virtual detector array, standard fihered backprojection, 
algebraic reconstruction techniques, and other tomographic imaging algorithms may be 
applied to generate image data of an object. 

In the examples shown here, the real detector array comprises three flat panel 
detectors arranged end-to-end, and angled to approximate an arc having a radius 

10 centered on the focal spot of the radiation source. It will be understood, however, that 
the principles of the invention can be used with actual detectors having any number of 
detector elements, including both ID line detectors and 2D panel detectors, where the 
geometry of the actual detector is neither equilinear or equiangular. In addition, the 
principles of the present invention can be advantageously employed in a system where 

1 5 one or more detectors are movable to various discrete positions along a line or arc 

relative to the x-ray source, such as described in co-pending U.S. Patent Application No. 
10/392,365, filed on March 18, 2003, the entire teachings of which are incorporated 
herein by reference. The principles of the present can also be used in a system in which 
the source and detector are tiltable about the focal spot of the source to obtain a larger 

20 field-of-view in the axial direction, such as described in co-pending U.S. application 
entitled "Cantilevered Gantry Positioning Apparatus for X-Ray Imaging System" 
(Attomey's Docket No. 3349.1004-001), filed on even date herewith, the entire 
teachings of which are incorporated herein by reference. 

While this invention has been particularly shown and described with references 

25 to preferred embodiments thereof, it will be understood by those skilled in the art that 
various changes in form and details may be made therein without departing firom the 
scope of the invention encompassed by the appended claims. 

Also, while the embodiments shown and described here relate in general to 
medical imaging, it will be understood that the invention may be used for numerous 
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other applications, including industrial applications, such as testing and analysis of 
materials, inspection of containers, and imaging of large objects. 



